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1.0 


INTRODUCTION  AND  SUMMARY 


This  is  a Final  Report  on  Contract  N00014-76-C-0721  entitled 
"IHR  Low  Altitude  Temperature  Profile  Retrieval."  The  program  is 
aimed  at  the  development  of  radiometric  techniques  for  the  remote 
measurement  of  atmospheric  vertical  and  slant  path  temperature 
and  humidity  distribution.  The  retrieved  profiles  are  expected 
to  be  instrumental  in  decreasing  the  errors  in  measured  target 
position  resulting  from  refraction  index  gradients  along  an 
atmospheric  propagation  path. 

Remote  profile  measurements  of  temperature  and  atmospheric 
gaseous  constituents  have  been  previously  implemented  by  radio- 

j 

metric  m.easurement  of  the  upwelllng  radiance  using  instruments 
mounted  on  aircraft,  satellite  and  balloon  platforms  (References 
1 and  2).  Unfortunately,  altitude  profile  resolution  obtainable 
with  similar  techniques  when  applied  to  up-looking  sensors  is 
inherently  poor  due  to  erratic  (non-peaked)  responses  of  the 
kernel  (weighting)  functions  associated  with  downwelling  radia- 
tion. The  information  thus  obtained  lacks  sufficient  detail  to 
retrieve  sharp  gradients  of  temperature  and  humidity  in  the  lower 
atmosphere  which  grossly  affect  local  variations  of  refractive 
index. 

An  alternate  technique  for  obtaining  more  accurate  tem.per- 

! 

ature  and  humidity  profiles  with  a grazing-angle  up-looking  radiom-  ' 

eter  has  been  investigated  and  is  discussed  in  detail  in  this  report. 

The  profiling  technique  is  based  on  measuring  radiation  from  pairs  i 

of  rotational  transition  absorption  lines  characteristic  of  abundant 


gases  (such  as  CO^ ) in  the  atmosphere.  The  technique  enables  the 
generation  of  difference- kernel  functions  which  are  well  suited  to 
the  task.  Spatially  selective  information  thus  generated  may  be 
utilized  for  high  accuracy  retrievals  of  temperature  and  humidity 
gradients  at  the  lower  altitude  atmospheric  levels  such  as  within 
the  Marine  Boundary  Layer. 

Retrievals  of  a double-inverted  temperature  profile  using 
the  differential  kernel  method  in  the  lower  atmosphere  have  been 
simulated  on  a computer.  For  this  effort  the  IJ  qm  CO^  band  is 
used.  The  model  profile  consisted  of  a double  inversion  at  575 
meters  and  at  75O  meters  above  the  ocean  surface.  An  iterative 
method  for  solutions  of  the  radiative  transfer  equations,  requir- 
ing no  apriori  statistical  information  on  atmospheric  parameters 
was  used  for  profile  retrieval.  The  general  methodology  of  the 
simulated  retrieval  process  is  discussed  in  the  report. 

A profile  retrieval  based  on  three  information  channels 
was  first  simulated.  The  retrieved  profile  obtained,  clearly 
reproduces  the  two  inversions,  however,  the  fidelity  of  the 
temperature  gradient  is  somewhat  poor  at  the  lower  altitude 
levels . 

A second  profile  retrieval  based  on  four  information  chan- 
nels, with  higher  spatial  selectivity  kernels  near  ground  level 
was  then  implemented.  The  resulting  retrieval  is  in  good  agree- 
ment with  the  actual  temperature  distribution  along  the  assumed 
slant  path. 
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Two  spectro-radiometric  techniques,  suitable  for  high 
spectral  resolution  remote  sensing  applications,  were  investi- 
gated on  the  program;  namely.  Infrared  Heterodyne  Radiometry 
(IHR)  and  Fabry-Perot  Etalon  (FPE)  spectroscopy.  A sensitivity 
analysis  has  been  performed  so  as  to  assess  the  performance 
limits  achievable  with  each  radiometric  type.  The  results  of 
the  analysis  show  the  FPE  to  be  superior  in  performance,  prim- 
arily due  to  its  higher  optical  throughput.  As  shown  in  this 
report,  reasonably  high  signal-to-noise  ratios  (SNR)  are  poten- 
tially achievable  with  an  FPE  system.  Basic  spectroscopic  tech- 
niques for  achieving  high  spectral  resolution  using  FPE  config- 
urations are  briefly  discussed  in  the  report. 
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GENERAL  SCOPE 


In  general,  variations  in  the  atmospheric  refraction  index 
may  be  classified  as: 

« Small  scale 
• Large  scale 

Small  scale  variations  are  associated  with  clear  air  turbu- 
lence phenomena  in  the  lower  atmosphere  and  are  manifested  in  random 
fluctuations  of  amplitude  and  phase  of  energy  transferred  through 
the  medium,  with  time  constants  in  the  millisecond  range.  The 
large  scale  refraction  index  profile;  of  interest  in  this  study, 
is  a slowly  varying  property  of  an  atmospheric  path  and  a function 
of  average  values  of  local  pressure,  temperature  and  humidity  along 
the  path. 

The  refractive  index,  n,  of  the  atmosphere  is  conveniently 
described  in  terms  of  the  refractivlty  N,  in  parts  per  million  by 
the  expression: 

n = 1 + 10"^  N (1) 

The  refractivlty  at  optical  wavelengths  is  related  to  local 
temperature  and  pressure  by: 

N::a77.6  ^ + 3.73  X 10^  ^ (2) 

where:  P = ambient  pressure  in  millibars 

T = temperature  in  degrees  Kelvin 
e = water  vapor  partial  pressure  in  millibars 
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The  downward  bending  curvature  of  a ray  propagating  at  a 
given  altitude  level  is  given  by  the  local  gradient  of  refractivity : 


77.6  <?P  , 

~T-  -^h  + 


+ 7.45  X 10 


Pt  , 3.73  X 10^  9e 

<?h 


(3) 


The  effect  of  humidity  is  appreciable  for  RF  and  microwave 
propagation.  At  optical  frequencies  however,  the  humidity  effect 
is  negligible  and  only  temperature  variations  must  be  considered. 

Absolute  values  of  local  temperature  humidity  and  pressure 
as  a function  of  altitude  entering  in  equation  (3),  may  be  estimated 
within  a small  error  from  reported  data  for  given  geographic  area 
and  measured  ground  (or  ocean)  temperature.  Sim.ilarly,  close  estim- 
ates on  pressure  gradients  are  obtained  from  universal  curves.  The 
situation  however  is  quite  different  for  the  case  of  local  tempera- 
ture and  humidity  gradients.  At  the  lower  altitudes  changes  in  the 
gradient  direction  (i.e.,  inversions)  of  temperature  occur  frequently 
as  a result  of  m.any  variables  (e.g.,  position  of  the  sun,  overcast 
condition,  surface  and  ambient  tem.perature,  wind  velocity  and  direc- 
tion). Likewise,  frequent  variations  occur  in  the  gradients  of  the 
humidity  profile  at  the  lower  atmosphere. 

Measured  temperature  gradients  at  lower  altitudes  using 
radiosonde  techniques  are  shown  in  Figure  1.  The  double  inversion 
in  the  gradient  ( superadlabatlc ) is  typical  for  mid-day  condition 
whereby  the  surface  temperature  is  higher  than  the  surrounding 
ambient  temperature.  The  neutral  gradient  shown  in  Figure  1 is 
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typical  for  early  morning  and  later  afternoon  conditions.  In  addi- 
tion a positive  inversion  at  ground  level  may  exist  at  night  when 
the  earth  is  cooling.  It  is  clear  that  the  uncertainty  in  the 
changing  gradients  of  temperature  in  the  vicinity  of  ground  or 
ocean  surfaces  along  with  changing  humidity  gradients  are  the 
largest  potential  sources  of  error  in  estimating  the  line  of  sight 
refractive  index. 

The  purpose  of  this  study  is  to  develop  techniques  for 
determining  the  fine  structure  temperature  and  humidity  gradient 
profiles  at  the  lower  altitudes  by  means  of  high  spectral  resolu- 
tion radiometric  measurements.  The  profile  information  thus  ob- 
tained should  make  it  feasible  to  assess  target  error  positions 
in  real  time  systems. 
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3.0  BASIC  PHENOMENOLOGY  IN  RADIOMETRIC  PROFILING 


The  average  apparent  radiance  sensed  by  a radiometer  viewing 
the  atmosphere  is  governed  by  the  steady-state  atmospheric  radiation 
transfer  equation  as  illustrated  schematically  in  Figure  2. 

For  a downlooking  radiometer  (Figure  2A),  the  sensible 
radiance  (neglecting  scattering)  is  given  by: 


o 


where : 


B , = Planck's  radiance  at  ground  and  i^^  altitude 
level,  respectively 

T , = atmospheric  transmission  coefficients  in  the 
total  propagation  path  from  ground  level  to 
receiver,  and  i^^  altitude  level  to  receiver, 
respectively 


9^i 

TfT 


= kernel  functions  of  the  transfer  equations 


'The  first  term  in  the  right-hand  side  of  equation  (4)  is  due  to  the 
earth  (or  ocean)  radiance  as  attenuated  in  its  total  propagation 
path.  The  second  term  represents  re-emission  of  radiance  in  the 
various  atmospheric  layers  within  the  propagation  path. 

The  transfer  equation  for  an  up-looking  radiometer  (Figure 
2B),  is  given  by  the  integral: 


^A  = 


-f  B,(T) 


dh 


(5) 


00 


where  the  radiance  is  due  to  emission  from  the  total  atmosphere, 
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The  method  of  fine  structure  profiling  consists  of  imple- 
menting high  spectral  resolution  sensor-channels,  ea’ch  selectively 
sensing  emitted  radiation  from  the  constituent  gas  emission  lines. 

As  shown  in  Reference  3,  the  spatial  resolution  of  a retrieved 
profile  based  on  radiometric  measurements  is  of  the  order  of  the 
spatial  width  of  the  kernel  function  in  the  spectral  sensor  channels. 

An  ideal  up-looking  radiometric  system  is  depicted  in  Fig- 
ure 3.  As  shown,  it  consists  of  a multitude  of  sensor  channels, 
each  having  idealized  rectangular  narrow  response  characteristics, 
sensitive  to  radiation  from  discrete  altitude  le  \e1s . These  ideal- 
ized weighting  functions  would  cause  each  channel  to  respond  only 
to  the  radiation  emitted  in  its  respective  altitude  level,  h. . 

The  actual  channel  responses  however,  are  normally  broadened  by 
natural  characteristics  of  the  gaseous  emission.  Hence,  the 
unique  altitude  behavior  is  not  available  and  each  channel  indeed 
senses  radiance  from  all  altitudes,  A crucial  result  of  this  effort 
therefore,  is  to  obtain  well  behaved,  peaked,  weighting  functions. 

Basic  parameters  associated  with  radiometric  profiling  are 
defined  in  Figure  4,  where  a down-looking  radiometer  senses  the 
radiation  from  an  "atmospheric  layer"  within  its  field-of- view. 

The  pertinent  radiometric  parameters  are  as  shown: 


B^(T)  = 
w^(p)  = 
€'(w,p)  = 
T = 
P = 


Planck's  radiance 

partial  pressure  of  active  gas 

gas  emissivity 

temperature 

ambient  pressure 
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Figure  3*  Ideal  Profiling  Rad"'  ometer  Concept 


r 

It  should  be  noted  that  though  gas  partial  pressure  and  emissivity 
are  also  functions  of  temperature,  this  functional  dependency  is 
small  as  compared  to  that  due  to  ambient  pressure,  and  is  neglected 
in  the  present  discussion. 

In  the  following  discussion  we  examine  the  difference  in 
shape  characteristics  associated  with  weighting  functions  for  up- 
looking  and  down-looking  radiometers. 

th 

The  apparent  radiance  emanating  from  the  i layer  as  sensed 
by  the  radiometer  shown  in  Figure  4 is  given  by: 

N.  = B.(T)  €^t.  (6) 

where  is  the  transmission  coefficient  of  the  total  path  from  i^^ 
level  to  receiver. 

To  relate  the  apparent  emission  expression  'to  the 

weighting  function  expression  we  resort  to  Kirchoff's  principle 
on  emission-absorption  equivalency.  Referring  to  Figure  4,  assume 
a radiance  source  Bg(T)  alternately  located  at  points  A and  B; 
i.e.,  below  and  above  the  atmospheric  layer  respectively.  The 
apnarent  absorption  within  the  layer  is  then  given  by: 


From  Kirchoff's  law  of  absorption-emission  equivalency,  the  multi- 
plying coefficients  of  Planck's  radiance  in  equation  (6)  is  equiva- 
lent to  the  transmissivity  difference  in  equation  (7),  which  may  be 

expressed  in  a differential  form,  namely: 

3x, 

A T - ^2  ~ — ^17  ® ) 
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are  the  kernels  of  weighting  functions 


The  spatial  derivatives 
and  are  linear  coefficients  of  the  temperature  dependent  Planck's 
radiance  in  equations  (4)  and  (5)«  spatial  selectivity  of  a 

radiometer  is  thus  determined  from  the  response  characteristics  of 


the  kernels  as  a function  of  altitude. 

The  difference  in  the  nature  of  the  kernel  responses  for  up- 
looking  and  down-looking  radiometers,  sensing  down  and  upwelling 
radiation  respectively,  may  be  visualized  by  referring  to  Figure  5. 
For  a down- looking  radiometer  (Figure  5A),  the  concentration  of  the 
active  gas  and  hence  its  emissivity  is  decreasing  in  consecutive 
atmospheric  layers  along  the  direction  of  radiation  (as  the  observer 
moves  closer  to  the  radiometer).  Ihe  transmission  coefficient  from 
each  consecutive  layer  on  the  other  hand,  increases.  Hence,  peak- 
ing of  the  product  €t  and  its  equivalent  weighting  function  will 
occur  somewhere  along  the  propagation  path.  For  an  up-looking 
radiometer  however,  no  peaking  will  occur  since  both  the  concen- 
tration and  the  transmission  increase  along  the  direction  of  radi- 


ation as  shown  in  Figure  5B. 

T5q)ical  altitude  dependence  of  monochromatic  weighting  func- 
tions are  shown  in  Figure  6 for  a down-looking  system,  and  Figure  7 
for  an  up-looking  system.  As  may  be  seen  from  the  figures,  unlike 
the  down-looking  case,  no  peaking  in  the  kernel  response  result  in 
up-looking  radiometers.  Hence,  the  spatial  resolution  of  a direct 
sensing  up-looking  radiometer  is  inherently  poor. 
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DIRECTION  OF  RADIATION 
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Figure  5.  Behavior  of  Kernels  for  Up-1  oking  and  Down-looking  Systems 
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DIFFERENTIAL  KERNEL  APPROACH 
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From  the  above  discussion  it  is  clear  that  the  weighting 
functions  resulting  in  spectral  channels  of  an  up-looking  radiom- 
eter based  on  direct  sensing  of  downwelling  radiation  suffer  from 
inherent  poor  spatial  selectivity.  This  suggests  that  additional 
properties  of  the  gaseous  emission  signature  must  be  utilized  in 
generating  selective  kernel  functions  for  up-looking  radiometers. 

The  method  developed  on  this  program  utilizes  the  steep  spectral 
slope  characteristics  of  single  rotational  line  signatures.  The 
technique  requires  sensing  radiation  from  sections  of  single  emis- 
sion lines  by  means  of  high  resolution  spectro-radiometers . The 
information  thus  obtained  enables  the  measurement  of  differences 
of  radiances  sensed  from  judiciously  selected  spectral  positions 
with  respect  to  the  centers  of  line-pairs  with  predetermined 
signature  strengths. 

The  basic  mathematical  i elationships  underlying  the  m.ethod 
of  generating  the  spatially  selective  difference  kernels  are  dis- 
cussed in  this  section.  The  general  expression  of  optical  depth 
as  a function  of  altitude  is  given  in  paragraph  ^.1,  and  from  it, 
m.athematical  expressions  of  kernel  functions  are  derived  in  para- 
graphs 4.2  and  4.3.  The  general  method  for  generating  the  spatially 
selective  functions  are  discussed  in  paragraph  4.4.  Plots  of  gen- 
erated spatially  selective  difference  kernels  for  carbon  dioxide 
and  water  vapor  are  shown. 


L 
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4.1 


ALTITUDE  FUNCTION  OF  A MONOCHROMATIC  ATMOSPHERIC  TRANSMISSION 
COEFFICIENT 


Figure  8 illustrates  a rotational  transition  line  shape,  with 
Lorentzian  attenuation  profile  representing  pressure  regions  in  the 
near  ground  levels.  The  monochromatic  absorption  coefficient  as  a 
function  of  frequency  is  given  by: 


«(v)  = ^ 


TT 


(v-v  Y 

^ O ' 


(9) 


where:  S = integrated  line  intensity 

*/■=  half-width  at  half  maximum  intensity 

v,v  = variable  and  line  center  frequency  (in  wavenumbers), 
respectively 

The  altitude  variation  of  the  absorption  coefficient  is  due  to  a 
linear  dependency  of  linewidth  with  pressure: 

'TqV  (10) 

where : 

p 

p = = normalized  pressure 

P, = pressure  at  altitude  and  ground  levels,  respectively 
= linev/idth  at  ground  level 

The  total  optical  depth,  Q(H)  from  altitude  H,  to  a ground 
based  radiometer  is  a function  of  gas  constituent  distribution  with 
altitude  as  well  as  variation  of  its  absorption  coefficient.  Thus: 


Q(H) 


a(h)  secG  w(h)  dh  = 


p(h)  w(h) 
C^+p^(h) 


dh 


(11) 


where : 

w(h)  = gas  content  per  unit  path  length  as  a function 
of  altitude 

v-v 

_ o = frequency  deviation  from  l^ne  center  normalized 
to  line  half-width  at  STP 
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Figure  8.  Lorentz  Line  Shape  Profile 


0 = zenith  angle  of  FOV 

p(h)  = normalized  pressure  as  a function  of  altitude 
The  transmission  coefficient  altitude  dependency  is  given  by: 

t(H)  = (12) 

4,2  MONOCHROMATIC  WEIGHTING  FUNCTION  RELATIONSHIPS  FOR  UNIFORMLY 
MIXED  GAS  SIGNATURES 

Signature  bands  of  uniformly  mixed  gases,  in  particular  the 

bands  of  carbon  dioxide  at  infrared  frequencies  are  utilized  for 

temperature  retrievals.  The  partial  pressure  of  a uniformly  mAxed 

gas  is  linear  with  pressure,  given  by: 


w(p)  - w^p 


(15) 


where:  w^  = absorber  content  at  STP  in  ( atm- cm )/k.m , The  variation 
of  pressure  with  altitude  is  expressible  by  an  exponential  relation- 


ship : 


(14) 


where:  h^  = scale  height  (^^.7,5  km) 
y = scale  height  parameter 

Substituting  equation  (I3)  into  equation  (11),  and  changing 
the  variable  of  integration  (from  equation  (14));  nam.ely: 


dp  = ^ e“  ^ dh  = - g-  dh  (I5) 

00  o 

the  optical  depth  of  equation  (11)  is  expressible  as  a function  of 
local  pressure,  namely: 
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where ; 


S w 


B = 


secO 


^ Vo 

p’  is  a variable  of  integration 

The  transmission  coefficient  is  now  readily  obtained  by  substituting 
the  expression  for  optical  depth  (equation  (l6))  into  equation  (12): 

B/2 


T = exp  J 


/ 

"^2  2 

B/2' 

r.2  2 

. In 

C +D 

C +p 

2 

1 + C^ 

P 

14C"^ 

^ “J 

> 

• ^ - 

(17) 


It  is  convenient  to  express  the  kernel  as  a spatial  derivative  with 
respect  to  the  scale  height  parameter  y;  thus,  from,  equations  (1?) 
and  (14), the  kernel  function  is  given  by: 


9p 

3P 


(1+C") 


(18) 


4.3 


MONOCHROMATIC  WEIGHTING  FUNCTION  RELATIONSHIPS  FOR  WATER 
VAPOR  SIGNATURE 


Water  vapor  distribution  in  the  atmosphere  varies  as  a func- 
tion of  humidity  conditions,  and  since  its  functional  relationship 
with  altitude  is  no  longer  linear,  the  transmission  coefficient  miay 
not  be  expressed  explicitly.  The  Integral  expression  of  the  trans- 
mission coefficient  as  a function  of  the  scale  height  parameter  is 
generally  given  from  equations  (11),  (12),  and  (14): 


exp  (-^y*  p(y)  dy) 


(19) 


whe  re : 


B = 


S h 


^ Vo 
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F(y)  = 


e-y 

+ e-2y 


w(y) 


w(y)  = HgO  content  per  unit  path  length  as  a function  of 
scale  height  parameter 


The  weighting  function  at  altitude  level  y^  is  then  given  by: 


3 T 

py 


3 F(yj_)  exp 


F(y)  dy 


(20) 


^ ^ GENERATION  OF  DIFFERENCE  KERNEL  INFORMATION  CHANNELS 


A typical  arrangement  of  channel  positions  for  generating 
difference  kernels  is  snown  in  Figure  9. 

To  first  illustrate  the  existence  of  a peaked  response  in 
difference  radiance  kernels,  assume  one  channel  is  tuned  to  line 
center  with  respect  to  line  1,  and  the  other  channel  at  the  far 
wings  with  respect  to  line  2.  Assuming  a uniformly  mixed  gas, 
referring  to  equation  (l8),  we  have; 


A.  For  the  first  channel; 


0=0 


— = 3,P^1 
py  1^ 


where  3,  = 


(21) 

absorption  coeffic- 
ient of  line  1 


C»l>  p 


(22) 


where  £*.re  the  absorption 

coefficient  and  normalized 
incremental  frequency  at 
line  2,  respectively 


The  difference  kernel  A in  the  two  channels  is 

A - R „2 

A — P-j  pi—  o P 


(23) 
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Figure  9.  Schematic  2- Channel  Arrangement  for  Generating  Difference 

Kernels 
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The  peak  of  the  difference  Kernel  is  obtained  by  minimizing 
the  function  (23);  namely: 


K 


2 p(ei-2) 


2 0. 


i- 


0 


(24) 


A minimum  of  the  A- functions  occur  at  infinity  (p  = 0),  and  a peak 
occurs  at  a pressure  level: 


■peak 


2 P, 


S 2c  2 

Pi  ^2 


l/(P,-2) 


(25) 


Thus  a peaked  difference  kernel  results  for  the  prescribed  ratios 
of  absorption  coefficients  and  channel  positions  in  the  two  lines. 

It  should  be  pointed  out  however,  that  the  difference  kernel  obtained 
from  channels  at  extreme  spectral  positions  are  rather  broad  due  to 
the  shallow  gradient  at  the  far  wing  of  the  line.  More  selective 
difference  kernels  result  for  a channel  position  at  a spectral 
region  of  steep  slope  in  the  absorption  coefficient  where  the 
sensitivity  to  pressure  variation  is  maximum. 

The  CO2  emission  band  at  I3  was  considered  for  generating 
differential  kernel  functions  for  temperature  profile  retrieval. 

The  technique  hov/ever  can  be  adapted  to  other  spectral  bands  at 
infrared  or  millimeter  wavelengths  where  discrete  signature  spectra 
due  to  uniformly  mixed  atmospheric  gases  exist. 

The  difference  in  kernels  for  two  channels  associated  with 
two  discrete  lines,  is  given  from  equation  (17)  by: 


25. 


1 


r 


1 

i 


i 


A - (^)i  - 


c^+v^ 

. i+c^^  J 


6^/2 


2 2 

C, 


-3. 


2._2 


Co^+P 


1+C. 


^2  2 
Co  +p^ 


(26) 


where:  ^’slate  to  signature  strength  of  the  two  lines 

^ as  defined  in  equation  (14) 

^1’  *"2  normalized  spectral  positions  of  the  two 
line  channels 


Line  signature  parameters  and  channel  positions  corresponding  to 
weighting  functions  with  peaks  at  prescribed  altitude  levels  were 
found  from  solutions  of  minimized  transcendental  equations  of  the 
form.: 


y 


y,- 


0 


(27) 


The  computations  were  performed  for  various  combinations  of  param- 
eters so  as  to  obtain  optimum  selective  and  contiguous  difference 
kernel  channels  for  altitude  levels  to  about  5 km.. 

Figures  10  and  11  show  difference-weighting  functions  obtained 
for  vertical  and  JO-degree  grazing-angle  up-looking  systems.  As 
shown,  .narrower  width  difference  kernels  result  at  regions  immedi- 
ately above  ground  level  and  broadening  toward  higher  altitudes. 

The  dotted  lines  shown  in  the  figures  represent  reversal  in  signs 
(foldovers)  which  will  tend  to  degrade  the  retrieval  accuracy.  In 
general,  appreciable  negative  response  values  in  the  vicinity  of 
ground  levels  are  undesirable.  To  that  end  the  kernel  functions 
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should  vanish  at  ground  level  which  is  accomplished  by  the  proper 
selection  of  the  relative  spectral  position  and  of  the 
channels.  Thus: 


A = 0 at  p = 1 

which  results  in  a constraint  on  the  absorption  coefficient  ratio 
given  by: 


1 + Cn 


1 + c. 


(28) 


Difference  kernel  functions  were  also  generated  for  water 
vapor  signature  lines  from  the  transmission  integral  equation  (20), 
which  are  given  by: 


0 1 r '^i  r 

1 ' "3y|2  " ^l^l^^i) 

o o' 

(29) 

T^ynical  humidity  distributions  for  standard  atmospheres  were  assumed 
(Reference  4),  and  corresponding  optical  depths  and  channel  positions 
were  computed  by  using  the  6-point  Gaussian  integration  method  on  a 
computer. 

Figure  12  shows  contiguous  peaked  difference-weighting  func- 
tions obtained  for  a 30-degree  grazing-angle  up-looking  radiometer, 
for  H^O  signature  lines  in  the  vicinity  of  4.7  |i,m. 

Results  of  computed  simulated  temperature  retrievals  using 
the  present  technique  are  discussed  below  in  Section  5. 
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Figure  12. 

DIFFERENTIAL  KERNEL  FUNCTIONS  FOR  WATER  VAPOR 
ALTITUDE  PROFILE  RETRIEVAL 


A 


5.0  SIMULATED  TEMPERATURE  PROFILE  RETRIEVALS 

Retrievals  of  a typical  superadiabatic  lapse  temperature 
profile  were  simulated  on  a computer  using  three  and  four  differ- 
ence-kernel information  channels.  The  basic  method  and  results 
obtained  are  discussed  in  the  following  paragraphs. 

5.1  RETRIEVAL  METHOD 

The  retrieval  algorithm  adapted  is  based  on  a modified  inverse 
solution  of  the  radiative  transfer  equation  in  conjunction  with  a 
relaxation  iterative  method  of  convergence  (References  5 and  6). 

The  advantage  of  this  method  over  the  inverse  m*atrix  retrieval 
techniques  is  that: 

a)  No  apriori  statistical  knowledge  on  atmospheric  param- 
eters are  required.  The  method  is  particularly  suited 
for  real  time  continuous  profile  updating  systems  at 
all  geographic  locations . 

b)  The  resultant  profile  is  generated  from  numerous 
solutions  at  small  Incremental  and  contiguous  alti- 
tude layers. 

The  basic  computer  routines  for  a temperature  simulated  retrieval 
analysis,  using  the  iterative  method,  are  shown  in  the  block  diagram 
of  Figure  13.  A somewhat  similar  computational  method  may  be  adapted 
for  humidity  retrieval. 

Referring  to  Figure  13,  the  basic  computational  steps  are  as 
follows : 


Figure  13. 

BLOCK  DIAGRAM  OF  COMPUTER  SIMULATED  PROFILE  RETRIVAL 


STEP  1.--A  temperature  altitude  profile  is  first  assumed 


and  designated  as  the  "true"  model  atmosphere. 

STEP  2. --The  radiance  differences  at  each  signature  channel 
corresponding  to  the  "true"  profile  are  computed  from  a digitized 
atmospheric  model  at  discrete  isobaric  incremental  regions,  and  are 
designated  as  "measured"  values. 

STEP  3. --An  arbitrary  temperature  profile  designated  as  a 
"guessed"  model  is  then  assumed,  serving  as  a starting  point  in 
the  iterative  computational  process.  In  consecutive  iterations 
the  model  is  updated  by  new  generated  profiles  in  Block  9. 

STEP  4. --The  "guessed"  (or  updated)  profile  is  digitized  at 
discrete  isobaric  layers,  and  difference  radiances  corresponding  to 
the  model  in  all  signature  channels  are  computed,  and  are  designated 
as  "comiputed"  values. 

STEP  5- --Difference  kernel  functions  at  each  discrete  spatial 
layer,  and  for  each  signature  channel  are  computed  based  on  the  known 
concentration  of  the  absorbing  gas  at  the  local  pressure  levels. 

STEP  6.--A  relaxation  constant  proportional  to  the  ratio  of 
"measured"  to  "computed"  values  is  set  for  each  signature  channel. 

This  constant  is  used  for  updating  the  model  profile  so  as  to  con- 
verge to  the  "true"  profile  in  a consecutive  step. 

STEP  7.--A  criterion  test  is  performed,  based  on  the  ratio 
value  in  step  6.  If  the  ratio  R in  all  signature  channels  is  close 
to  unity  the  iteration  process  is  terminated.  Otherwise,  the  process 
is  continued  at  step  8. 

i 

i 
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STEP  8. --The  discrete  temperatures  of  the  profile  model 
(step  4)  are  updated  in  proportion  to  the  relaxation  constant 
corresponding  to  the  given  signature  channel  (step  6).  In  essence 
new  updated  temperature  profiles  are  generated  at  each  signature 
channel . 

An  iterated  profile  is  finally  obtained  from  the  individual 
channel  models  (step  6),  which  is  based  on  an  average  weighting 
in  proportion  to  the  kernel  values  corresponding  to  the  channel 
at  each  altitude  level. 

The  program  has  been  written  and  executed  for  three  and  four 
difference  kernel  retrieval  simulations.  Tabulated  spectroscopic 
signature  parameters  from  Reference  7 were  used.  The  profiles  were 
retrieved  in  intervals  of  37.5  meters  for  the  first  75O  meters  and 
in  75  meter  intervals  from  75O  to  2.25  km  above  ground  level,  and 
750  meter  Intervals  above  2.25  km.  Between  30  to  50  iterations 
were  required  for  convergence. 

5.2  SIMULATED  RETRIEVAL  RESULTS 

Simulated  retrievals  for  an  assumed  superadiabatic  tempera- 
ture profile  with  a double  inversion  at  375  meters  and  75O  meters 
above  ground  level,  were  programmed  on  a computer. 

Figure  14  shows  a retrieval  based  on  a three  information 
channel  system,  with  weighting  functions  as  shown  in  Figure  11. 
Curve  A in  the  figure  is  the  assumed  "true"  atmospheric  profile, 
curve  B is  the  "arbitrary  assumed"  initial  model  in  the  iteration, 
and  curve  C is  the  retrieved  model  profile.  As  shown,  the  two 
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COMPUTER  SIMULATED  RETRIEVAL 
(3  INFORMATION  CHANNELS) 


inversions  are  reproduced  in  the  retrieved  profile  with  a dis- 
placement in  gradient  of  O.5  to  2^K  on  the  temperature  scale  at 
the  upper  altitude  levels  and  with  a somewhat  poorer  fidelity 
immediately  above  ground  level. 

A retrieved  profile  obtained  on  the  computer  using  four 
information  channels  with  spatially  narrower  width  weighting 
functions  for  channels  in  the  vicinity  of  ground  level,  is  shown 
in  Figure  Ip.  The  fidelity  of  the  profile  gradient  is  consider- 
ably improved  at  the  lower  regions. 

The  discrepancies  in  the  fine  structure  retrievals  are 
the  result  of  the  negative  foldovers  in  the  difference  kernels 
at  upper  altitudes.  These  introduce  errors  due  to  uncertain- 
ties in  "measured"  radiance  differences  associated  with  the 
information  channels;  an  effect  somewhat  similar  to  that  pro- 
duced by  system  noise.  Hence,  no  further  degradations  due  to 
noise  are  expected  to  occur  for  system.s  having  reasonable  high 
signal- to-noise  ratios. 
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COMPUTER  SIMULATED  TEMPERATURE  ALTITUDE 
PROFILE  RETRIEVAL  (4-CHANNEL  SYSTEM) 


6.0 


SYSTEM  SENSITIVITY  CONSIDERATIONS 


As  discussed  in  Section  5,  the  present  technique  of  selec- 
tive profiling  is  based  on  quasi-monochroraatic  measurements  of 
radiance  differences  at  prescribed  spectral  positions  with  respect 
to  absorption  line  pairs.  Thus  very  high  spectral  resolution  in 
the  radiometric  channels  is  required.  Hence,  there  are  two  key 
requirements  associated  with  the  spectro-radiometric  instrum.enta- 
tion,  namely: 

a)  High  spectral  resolution.  The  instrument  must  be 
capable  of  resolving  spectra  of  the  order  of  10”^ 
to  10"^  cm"^  (i.e.,  approximately  1 percent  of 
absorption  linewidth). 

b)  High  sensitivity.  The  instrument  must  be  capable 
of  sensing  small  signal  differences  from  narrow 
spectral  width  channels  with  sufficient  signal-to- 
noise  ratios. 

Two  candidate  systems  suitable  for  remote  detection  with  the  required 
high  resolution  were  investigated  on  the  program: 

a)  The  coherent  infrared  heterodyne  radiometer  (IHR) 

b)  The  Fabry- Perot  Etalon  (FPE)  spectro-radiometer 

In  the  IHR  technique,  a laser  local  oscillator  signal  (e.g., 
tunable  diode  laser),  is  coherently  mixed  with  an  incoming  source 
irradiance,  resulting  in  an  intermediate  frequency  (IF)  at  the  out- 
put of  the  photomixer.  Extremely  high  spectral  resolution  is  attain- 
able with  an  IHR  vfhich  ultimately  is  limited  only  by  the  stability  of 
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the  laser  source. 


r 

The  FPE  spectrometer  basically  consists  of  two  high  reflec- 
tive parallel  plate  mirrors,  whereby  high  resolution  bandpass 
response  characteristics  result  from  interferences  of  multiply 
reflected  beams.  The  spectral  resolution  in  the  FPE  is  ultimately 
limited  by  achievable  etalon  finesses  and  considerations  of  com- 
plexity in  practical  implementations  of  multi-etalon  structures 
required  for  rejecting  higher  order  passbands. 

An  analysis  on  relative  merit  and  performance  limits  of 
the  two  spectral  techniques  for  the  present  application  is  given 
below.  Key  spectral  techniques  associated  with  FPE  are  also 
discussed  in  this  section. 

6.1  IHR  SENSITIVITY 

The  sensitivity  of  a properly  designed  IHR  system  is  quantum- 
noise-limited  at  the  operating  optical  frequency.  The  noise  equiva- 
lent power  (NEP),  defined  as  the  input  signal  required  to  result  in 
a slgnal-to-nolse  (SNR)  of  unity,  is  hence  given  by: 

NEP  = hf  7—  ' (watts)  (30) 

I 

1 

, where:  f = optical  frequency  (Hz) 

i 2 

h = Planck's  constant  (watts-sec  ) 

Af  = IF  bandwidth  (Hz) 

T = video  integration  time  (sec) 

For  a spectral  signal  radiance  (AN),  resulting  from  a differ- 
ence signal  in  an  absorption  line  pair,  the  detected  signal  power  is 
given  by: 

Pp  = I ^ Ao  T^-n.  AN  AV  (31) 
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39. 


where : 


n = quantum  efficiency 

p ' 

= aperture  area  (cm  ) 

= optical  efficiency 

= solid  angle  field  of  view  (sr)  j 

P 1 i 

AN  = spectral  radiance  difference  (W/cm  -sr-cm~  ) j 

Av  = optical  bandwidth  (cm~^) 

and  the  factor  i accounts  for  polarization  selection  heterodyne 
rejection. 

For  optimum  efficiency  of  the  IHR,  the  throughput  defined  as  the 
product  A^,  is  diffraction  limited  requiring: 

n A^  = (52) 

where:  A = wavelength  (cm) 

From  relationships  in  equations  (3O),  (31)>  and  (52),  and  noting 
that : 

f = cv 

Af  = c (Av) 

. = 1 

''  V 

where:  c = speed  of  light  (cm/sec) 

an  expression  for  the  SNR  in  an  IHR  is  given  by: 


SNR  = 


D 

NEP 


AN  . 7] 


AV  T 


(33) 


2 hc^'  “ v 

We  define  a noise  equivalent  spectral  radiance,  ANj^^,  which  corres- 
ponds to  unity  SNR.  Thus,  from  equation  (33)»  setting  SNR  = 1 and 
solving  for  AN: 
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6.2  FPE  SPECTRO- RADIOMETER  SENSITIVITY 

The  sensitivity  of  a background  limited  (BLIP)  operated  FPE 
receiver  is  limited  by  fluctuations  of  background  photons  incident 
within  its  field-of-view,  and  within  its  spectral  acceptance  window. 
The  total  background  power  collected  by  the  FPE  radiometer,  viewing 
a background  scene  with  spectral  radiance  is  given  by: 

D 

Pfl  = " "o  •'o “b  " (’5) 

where : 

7?  = quantum  efficiency 

2 

Aq  = area  of  optical  telescope  (cm  ) 

Tq  = optical  efficiency 

il  = solid  angle  field  of  view  (sr) 

= background  spectral  radiance  as  a function  of 

P - 1 

temperature  and  emlssivity  (W/cm  -sr-cm  ) 

Av  = optical  bandwidth  (cm~^) 

The  corresponding  number  of  photons,  Q,  generated  in  the  detector 
over  the  duration  of  an  integration  time,  t,  is  given  by: 

Pb^ 

Q = (number  of  photons)  (3^) 

where:  hcv  = energy  per  photon  (W-sec) 

The  number  of  fluctuating  noise  photons,  q,  is  given  from  Poisson 
statistics  by  the  square  root  of  the  number  of  background  photons: 


q = ^ 


(57) 


and  the  corresponding  nol' e power  Integrated  over  time  t is  given 


= 2 hcv  q/t  (38) 

ivhere  the  factor  of  2 accounts  for  generation- recombination  noise 
processes , 

Combining  equations  (35)  through  (38)j  the  noise  power  is  given 
in  equation  (39): 


Pj^  = J^hcv  Ng  Av  7?  f 


(watts)  (39) 


For  a given  spectral  signal  radiance  AN,  resulting  from  a 
difference  signal  in  an  absorption  line  pair,  the  detected  signal 
•p'-v’er  is  given  by: 

Pg  = 77  Ap  t^JIAN  av  (to) 

Thus  the  SNR  is  given  from  the  above  equations  (39)  and  (tO): 

Pg  f 77  A T JIAV 

SNS  - ^ . an  I J (1.1) 

The  noise  equivalent  signal  spectral  radiance  ANj^g  is  obtained 
by  setting  SNR  = 1 and  solving  for  AN  in  equation  (41): 


4hcv  Ng 
T]  A^T^-n-Avt 


(watts  ) 


The  throughput  properties  of  the  FPE  will  now  be  considered. 
Since  Fabry-Perot  etalons  have  limited  fields-of-view  at  any  one 
resonant  condition,  energy  incident  appreciably  off  the  optical  axis 
will  propagate  with  difference  resonance  properties  in  the  etalon. 
The  limitation  on  FOV  0,  due  to  this  shift  is  given  by; 


0 = 2 (rad) 


Thus,  energy  incident  outside  the  angular  POV  0 will  not  contribute 
to  signal  power.  The  solid  angle  ^ in  equation  (42)  is  thus  con- 
strained by  © from  equation  (43),  and  since: 

-a  = J ©2  (sr)  (44) 

the  dependency  of  AN™  on  optical  bandwidth  is  given  by: 
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6.3  COMPARISON  OF  SENSITIVITY  LIMITS  IN  AN  IHR  AND  FPE  SPECTRO- 

RADIOMETER 

A comparison  of  ultimate  sensitivities  achievable  in  the  two 
spectro-radiometric  systems  is  showTi  in  Figure  l6.  The  values  plotted 
are  based  on  equations  (34)  and  (45)  for  the  IHR  and  FPE  respectively 
and  are  normalized  to  a radiance  of  300°K  blacEbody,  thus  represent- 
ing the  minimum  detectable  fractional  change  in  an  arbitrarily  assumed 
ambient  background  v/ith  unity  emissivity. 

As  may  be  seen  from  Figure  l6,  appreciably  greater  sensitivity 
results  with  a FPE  radiometer  as  compared  to  the  IHR,  which  is  par- 
ticularly pronounced  at  the  shorter  wavelengths  where  the  IHR  sensi- 
tivity is  inherently  lower.  The  inferior  performance  of  the  IHR 
results  primarily  due  to  the  small  throughput  associated  with  a 
coherent  receiver  which  requires  a diffraction  limited  FOV.  As  may 
be  noted,  the  minimum  detectable  spectral  radiance  values  plotted 
in  the  figure  are  normalized  to  the  square  root  of  the  optical  band- 
width. Since  the  throughput  in  an  FPE  radiometer  is  a function  of 
optical  bandwidth,  several  parametric  curves  are  shown  in  the  figure. 


figure  16.  COMPARATIVE  SENSITIVITY  OF 
SPECTRO-RADIOMETIC  SYSTEMS 


6.4 


ESTIMATED  PERFORMANCE 


An  estimated  sensitivity  level  of  performance  of  an  FPE 
soectro- radiometer  implemented  for  temperature  profiling  is  shown 
in  Figure  I7.  The  assumes  system  parameters  are: 

Diameter  of  optical  telescope  = 4" 

Optical  bandwidth  = 3 x 10~^  cm~^  (100  MHz) 

Quantiim/optical  efficiency  = 50  percent 
Observation  time  per  measurement  = 5 seconds 
The  bottom  curve  in  Figure  17  represents  the  minimum  resolvable 
fractional  change  in  ambient  (assumed  here:  T = 300°K,  1),  corres- 

ponding to  SNR  = 1.  The  shaded  area  represents  the  range  in  the  dif- 
ference signal  radiances  resulting  from  pairs  of  absorption  lines 
corresponding  to  various  spatial  selectivity  channels.  As  may  be 
seen,  reasonably  high  SNRs  are  obtainable  with  the  system. 

Figure  18  shows  required  D*  as  a function  of  operating  wave- 
length and  the  corresponding  number  of  photon  flux  densities.  Shovm 
detectivity  levels  are  within  the  state  of  the  art  of  currently 
available  detectors. 

6. 5 BASIC  FPE  SPECTROSCOPIC  INSTRUMENTATION 

A typical  response  of  an  FPE  is  shown  in  Figure  19.  As  shown, 
it  consists  of  a multiplicity  of  periodic  responses,  each  designated 
by  an  order  m^.  The  frequency  spacing  between  two  consecutive  band 
passes,  referred  to  as  the  free  spectral  range  is  a function 

of  the  optical  path  length  between  the  two  plates  as  shown  in  the 

equation  in  the  figure.  The  ratio  of  the  bandwidths  to  the  free 

✓ 

spectral  range  is  given  by  the  etalon  finesse  ^ , which  for  ideal 
etalons  is  a function  of  plate  reflectivity,  however  in  practice 
is  often  limited  by  the  achievable  plate  parallelism. 
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REQUIRED  D*  cm  Hz 


(PHOT/cm^  sec) 


m-1  m m+1  m+2 


FABRY  PEROT  ETALON (FPE) RESPONSE  CHARACTERISTICS 


A block  diagram  of  a basic  FPE  configuration  for  high 
spectral  resolution  measurements  is  shown  in  Figure  20.  As  shown, 
it  consists  of  a dielectric  interference  filter  followed  by  several 
etalon  filters  each  with  a decreasing  free-spectral  range.  Figure 
20b  shows  single  line  output  with  an  idealized  bandpass  response 
of  the  dielectric  filter.  To  that  end,  the  free  spectral  range 
(FSR)  of  the  first  etalon  is  made  equal  or  greater  than  the  width 
of  the  interference  filter,  so  as  to  assure  a single  order  bandpass 
response  within  the  filter  bandwidth  (see  Figure  20c).  The  ultim;ate 
high  spectral  resolution  is  obtained  from  consecutive  etalons,  each 
having  its  FSR  equal  to  the  bandwidth  of  the  preceding  etalon  as 
shown  in  Figures  20d  and  20e.  The  number  of  required  FPE  stages 
in  the  configuration  is  a function  of  the  achievable  finesse  for 
a desired  spectral  resolution. 
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FILTER  ^ 


PROPOSED  BASIC  SPECTROSCOPIC  INSTRUMENTATION 
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